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Abstract 

It can be expected that the respective endpoints of the Gregory-Laflamme black brane 
instability and the Rayleigh-Plateau membrane instability are related because the bifur- 
cation diagrams of the black hole-black string system and the liquid drop-liquid bridge 
system display many similarities. In this paper, we investigate the non-linear dynamics of 
the Rayleigh-Plateau instability in a range of dimensions, including the critical dimension 
at which the phase structure changes. We show that near the critical dimension and above, 
depending on a parameter in initial conditions an unstable cylinder will either pinch off 
or converge to an equilibrium state. The equilibrium state is apparently non-uniform but 
has a constant mean curvature everywhere. The results suggest that in the gravity side, 
near the critical dimension and above, the final state of an unstable black string (which is 
not too long) is a non-uniform black string. The equation of motion adopted to describe 
the dynamics is the surface diffusion equation, which was originally proposed to describe 
a grooving process of heated metal surfaces. An interesting correspondence between the 
diffusion dynamics and black hole (thermo)dynamics is discussed. 



1 Introduction 



A variety of black hole solutions exist in higher dimensional spacetimes even with the asymptotic 
flatness condition [1]. Furthermore, the compactification of an extra dimension has been known 
to make the phase structure richer [2,3]. For the spacetimes in which one spatial dimension is 
compactified on a circle, three phases of black objects are known, i.e., a localized black hole, 
a uniform black string, and a non-uniform black string [4-6]. Although the phase diagrams, 
which are now available in 5 and 6 dimensions in the best form [6], suggest that the transition 
between them would occur, our knowledge about stability/dynamics of each phase is quite 
restricted. What we know is that the uniform black strings are unstable perturbatively if the 
linear dimension is too large as shown by Gregory and Laflamme [7]. One might expect the 
endpoint of the instability to be a localized black hole since it has a larger entropy than the 
original uniform black string. Horowitz and Maeda, however, argued that the such a topology 
changing transition cannot occur due to the "no-tear" nature of the horizon and that the 
endpoint would be an non-uniform black string [8] (see also [9]). Then, a general relativistic 
simulation was followed [10,11] and some suggestions were obtained, but no one has seen the 
unstable black string to pinch off nor any symptom of the convergence to an equilibrium state. 

One of the most interesting aspects of the above black hole-black string system is that the 
phase structure seems to change drastically at a threshold spacetime dimension, d = : = 
14 [12]. In d < d*, the uniform black string is thought to transit suddenly to a localized black 
hole (or probably also to a non-uniform black string, depending on dimension) as the compact 
dimension is adiabatically stretched with the mass kept fixed [5, 13]. While in d > d*, the 
uniform black string is thought to transit smoothly to an infinitesimally deformed non-uniform 
black string, which would be stable unlike that in d < d*. Thus, it is quite interesting to 
understand phase structures above/around the critical dimensions and to know whether or not 
the stability of non-uniform black string indeed changes around the critical dimension. We note 
that the number of a critical dimension depends on the type of ensemble [14], momentum [15], 
and charge [16]. For example, the critical dimension in a canonical ensemble is d* = 13. 

The membrane paradigm of black holes [17] often provides useful tools to investigate how 
black holes actually behave and a new perspective for the black hole physics [18]. In the con- 
text of the black hole-black string system, a lot of similarities between the Gregory-Laflamme 
instability and the Rayleigh-Plateau instability, which is a universal instability of extended flu- 
ids/membranes [19,20], were pointed out perturbatively [21-23]. Then, the previous paper [24] 
by the present author and a collaborator revealed that the similarities of the phase structures 
between the black hole-black string system and a liquid drop-liquid bridge system persist up 
to non-linear regimes. In particular, it was shown that a critical dimension similar to that in 
the gravity side [12] exists in the fluid side. That is, the phase diagram implies that a Uniform 
Bridge (UB) suddenly transits to either a Non-Uniform Bridge (NUB, known as the Delaunay 
unduloid [25]) or a Spherical Drop (SD) below a critical space dimension, D < D* := 12, as the 
compact dimension is adiabatically stretched with the volume kept fixed. While in D > D*, 
the phase diagram implies that the UB smoothly transits to an infinitesimally deformed NUB 
just like the transition of the uniform black string to non-uniform black string in d > d*. 

In the previous paper [24] , the equilibria of a liquid were obtained by the variational calculus 
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which minimizes/extremalizes the surface area of the liquid while keeping the volume kept fixed 
(i.e., Plateau's problem or the capillary minimization problem). Thus, a liquid surface obtained 
is a minimal surface, which has a constant mean curvature everywhere. The next step would 
be to establish the (in)stability of each phase. Although one can investigate the local stability 
by the second variation of an area functional, it may be more interesting to see directly the 
global (in) stability adopting a suitable dynamical equation of motion. What kind of dynam- 
ics is appropriate? A necessary condition for a candidate is to have the equilibria represented 
by the constant mean curvature surfaces. In addition, it should be required that the motion 
decreases the surface area while keeping the volume fixed throughout time evolutions. Is there 
such a convenient dynamical equation? Yes, the surface diffusion equation, which was originally 
proposed to describe a grooving process of heated polycrystal surfaces in metallurgy [26], does 
have all properties raised above. Within the framework of the surface diffusion, it is known that 
in the usual 3 dimensional (3D) space the SD is stable, the NUB is unstable, and the UB is un- 
stable (the Rayleigh-Plateau instability) if its linear dimension is longer that its circumference. 
It is also known that the unstable UB dynamically transits to the SD via a topology changing 
transition in a self-similar manner [27, 28] (see [28] for a comprehensive study of axisymmetric 
dynamics) . 

In this paper, we investigate the linear and non-linear properties of the Rayleigh-Plateau 
instability and especially their endpoints in the surface diffusion dynamics. We will show (i) that 
the Rayleigh-Plateau instability has a dimensional dependence similar to that of the Gregory- 
Laflamme instability; (ii) an interesting aspect of the Rayleigh-Plateau instability appearing in 
higher-order perturbations; (iii) that an unstable UB pinches off in lower dimensional spaces 
(e.g., in 4D), while the unstable long (but not too long) UB converges to the NUB in higher 
dimensional spaces (e.g., in 12D); and point out (iv) a remarkable correspondence between the 
surface diffusion dynamics and the black hole thermodynamics. 



2 Axisymmetric Surface Diffusion 
2.1 Basic Equation 

The surface diffusion is a mass transport phenomenon on material surfaces driven by the gra- 
dient of a curvature [26] . A superficial flux of particle numbers is given via the Nernst-Einstein 
relation by 

J = -AV s k , (1) 

where A is a constant, V s the surface Laplacian, and k the mean curvature of the surface. Local 
volume conservation results in the motion of the surface: 

(m + V s • J = , (2) 
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Figure 1: An axisymmctric hypcrsurfacc in (n + 2)-dimcnsional space 



where u is the normal velocity of the surface and p is a constant density of atom numbers. 
Thus, the motion of surface obeys the surface diffusion equation^] 



u 



p~ X AV*K =: BA s k . (3) 



From now on we concentrate on the motion of axisymmetric (n + l)-dimensional "metal 
surfaces" in D := (n + 2)-dimensional space (n > 1). We denote the instantaneous radius of 
S n at z by r = r(t, z), where t is the time and z is the axis of symmetry (see Fig. [1]). In this 
coordinates, the mean curvature of the surface is given by 



n r" 



^VTT^ (1 + r' 2 f' 2 



(4) 



where X' = d z X. See Appendix for the calculation detail. The first and second terms in Eq. (J3j) 
have the meanings of the azimuthal and axial principal curvatures, respectively. The normal 
velocity and the surface Laplacian are given by 

- = ^SL= , A s = - ±== d z (^—d z ) . (5) 



r 



The combination of Eqs. ([3"lh pi), and (jSJ) gives the basic equation governing the dynamics of 
axisymmetric surface diffusions 



d t r = ^-d z I d z 



r n Vvl + r 



ii 



n r" 



(1 + r' 2 ) 3 / 2 J 



(6) 



which is first order and fourth order with respect to t and z, respectively. It is important 
to notice that from Eq. (J3j) the constant mean curvature surfaces (k = const.), e.g., the Uni- 
form Bridge (UB), Non- Uniform Bridge (NUB), and Spherical Drop (SD) obtained in [24], are 
stationary solutions of the surface diffusion equation. 



2.2 Volume Conservation and Area Decreasing 

Now, let us see two important aspects of the surface diffusion dynamics, i.e., the conservation 
of volume and the decreasing of surface area. When we consider an object extending in a region 

1 In the original theory of the thermal grooving [26], the constant B is given by B = Ds^i^is/ksT, where D s 
is the coefficient of surface diffusion, 7 the surface-free energy per unit area, v = p~ x the molecular volume, v 
the number of atoms per unit area, and T the temperature. 



of z G P := [-L/2, L/2], the instantaneous volume and surface area of the body can be written 

as 

V{t) = n n J dzr n+1 (t,z) , 

A(t) = (n + l)fi n J dzVl + r' 2 r n (t, z) . (7) 

Here, VL n := 7r (n+1 )/ 2 /T[(n + l)/2 + 1] is the volume of a unit n-sphere. Taking a time derivative 
of these equations, one can show 



V(t) = B(n + l)n r 



VT4 
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dP 
n,j2 



r n K 12 . r n (BKK' + r'd t r) 

A(t) = -B(n + l)n n / dz - +{n + l)Q n 



p 



.12 



(8) 

dP 



From equations in (jSj), we see that the volume does not change (V = 0) and the surface area 
decreases (A < 0) throughout a time evolution if the following periodic boundary conditions 
hold, 

r« (t, -L/2) = r« (t, L/2) , «« (t, -L/2) = «« (t, L/2) , (Z = 0, 1) . (9) 

Here, := The equality A = holds when the hypersurface has a homogeneous mean 
curvature, k' = 0. 

2.3 Rayleigh-Plateau Instability 

The linear perturbation around a UB sets a starting point of the investigation of dynamics. We 
expand r(t, z) around the UB as r(t, r) = ro + sri(t, z) with a small expansion parameter e and 
a radius of UB r . Substituting this expansion into Eq. (jHJ), we have the linear perturbation 
equation at 0(e), 



ii 



( + r i ] ) 



Now, we consider a superposition of N sinusoidal waves as an initial condition given by 



ri(0, z) = ^ cos(fcj2;) , (11) 



where aj and fc, are constants. Then, the perturbative equation fflOl) is solved to give 



iV 



r 1 (t,z) = ^a i e^)*cos(M , (12) 

i=l 

u(k):=- 4 (kr ) 2 [n-(kr ) 2 } . (13) 
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Figure 2: Growth rate of the Rayleigh-Plateau instability from 3D (bottom curve) to 6D (top curve). 



Equations (1T21) and ( {TBI tell us that the initial perturbation of wavenumber fc, grows exponen- 
tially if < ki < k CT := \fnjrQ for a given r , which implies the Rayleigh-Plateau instability 
is a long wavelength instability just like the Gregory-Laflamme instability. Furthermore, the 
dimensional dependence of this critical wavenumber is similar to that of the Gregory-Laflamme 
instability, and they indeed coincide each other in the large dimension limit, n — > +oo [29,30] 
(see also [21]). We plot the dispersion relation ( 113]) in Fig. [2j The peak of growth rate, 
w m£K := Lu(y^n/2/r ) = Bn 2 /Ar^, also increases with dimension, which is similar to the Gregory- 
Laflamme instability again [7]. One difference between two instabilities seems to exist near the 
infrared region kr$ -C 1, where the growth rate of the Rayleigh-Plateau instability is suppressed 
faster than the Gregory-Laflamme instability. This behavior would come from the "short dis- 
tance nature" of the curvature driven surface diffusion, relating to the fourth derivative term in 
Eq. (flOl) . Note that the dimensional dependence of the critical wavenumber, k CT oc \/n, comes 
not from the Laplacian in Eq. (jSJ) but from the azimuthal principal curvature in Eq. (j3J). There- 
fore, we can say that the dimensional dependence comes from that of the background curvature 
of S n rather than that of the "interaction" represented by the higher-order derivatives. 

One can generalize the above linear perturbation of UB to higher-order ones. In Appendix, 
we present the result of second-order perturbation. Here, we just mention an interesting aspect 
of Rayleigh-Plateau instability peculiar to the non-linear regime. If we set ki > y/n/r , (1 < 
% < N) in the initial condition (fTTj) . the linear perturbation analysis tells us that the initial 
perturbations will be exponentially suppressed. Even in this case, however, a certain class 
of modes in the second-order perturbation grows if there exists a couple of sinusoidal waves 
satisfying \k{ — kj\ < vW r o ; 7^ j) i n the initial perturbation [27]. 
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3 Non-Linear Dynamics 



3.1 Boundary and Initial Conditions 

Now, we are ready for the numerical investigation to know the endpoint of the Rayleigh-Plateau 
instability. The partial differential equation ([6]) is solved for various n. We solve the equation in 
the region of z G P = [-L/2, L/2] with the periodic boundary condition (Q, which is equivalent 
to 

r«(*, -L/2) = r®{t, L/2) , (1 = 0, 1, 2, 3) . (14) 

With boundary condition (fT4l) . one can model the motion of an infinitely long periodic body. 
For simplicity, a UB given a single sinusoidal perturbation is chosen as the initial condition: 

r(0,z) = r [1 + e cos (pz)} , p = 2n / L , L = (1 + 8)L cr , (15) 

where L cr := 2n/k cr is the critical wavelength, e and 5 are small numbers (5 > 0). Thus, 
we have a couple of free parameters (e, 5) in our minimal-setting simulation. In the following 
numerical calculation, we will set these parameters to 

£ = -0.10, 5 = 0.12. (16) 

The positivity of 5 implies that the initial perturbation has a wavenumber smaller than the 
critical wavenumber (p < k cr or equivalently L > L cr ) and will grow at least during an early 
time of the evolution. A typical dynamical time scale can be defined using the growth rate of 
linear perturbation ([TBI as 

t dyn := l/w(p) . (17) 
For example, one can normalize length and time by r and tdyn, respectively. 




Figure 3: Time evolution of r = r(t, z) in 4D (left) and 12D (right). 
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Figure 4: Snapshots at selected moments of time in 4D (top) and 12D (bottom). The evolution results in 
pinching off in 4D and convergence to a non-uniform equilibrium state in 12D. The final equilibrium (12D) is 
apparently highly non- uniform but has a constant mean curvature everywhere (see Fig. [5]). 



3.2 Fate of Instability in 4D and 12D 

We will show typical behaviors taking two examples, D = n + 2 = 4 (< D*) and D = n + 2 = 12 
(= -D*). The time evolution of r = r(t,z) obtained numerically is shown in Fig. [3] for these 
cases. 

In 4D, the given initial perturbation continues to grow and the numerical solution suggests 
the breakup of the body to become a pair of separated drops near t/t dyn ~ 0.92, just like 
the behavior observed in the 3D axisymmetric surface diffusion [27,28]. Since the breakup, 
represented by r = 0, is a singularity in our coordinates, we have to change variables before 
the breakup. Although the continuation of calculation through the breakup is possible with 
some prescriptions [27,28], we simply stopped the calculation since to know the detail of the 
breakup is not the purpose of this paper. Nevertheless, we can argue that the numerical result 
implies the pinch off from the behavior of mean curvature, as we will see soon. In 12D, the 
given initial perturbation grows during an early time, which is consistent with the linear- and 
second-order perturbations. However, the numerical solution suggests that the surface begins 
to converge to a certain non-uniform configuration due to non-linear effects. Several snapshots 
of the hypersurface at selected moments of time are shown in Fig. HI 

Now, we see the evolution of the mean curvature k. The time evolution of k is shown in 
Fig. [5] for both dimensions. Note that in both cases the initial sinusoidal perturbation breaks 
the constancy of the mean curvature of the background UB, Ktjb = n/r Q . In 4D, the mean 
curvature continues to grow with an increase of spatial gradient especially around the center, 
z = 0. This means that the current |«7| continues to grow and the surface deviates from an 
equilibrium state. This behavior supports our interpretation that the surface pinches off and 
also approaches the other equilibrium state, i.e., the SD which has a smaller area [24]. While in 
12D, the mean curvature becomes a global constant after a few oscillations localized near the 
center. This behavior implies the surface converges to the equilibrium state having a constant 
mean curvature, i.e., the NUB obtained in the previous paper [24]. 
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Figure 5: Time evolution of the mean curvature k = n{t,z) in 4D (left) and 12D (right). Note that the 
viewpoint is different from that in Fig. 



3.3 Initial Condition Dependence and Fate in Other Dimensions 

To state criteria for a body to pinch off or converge to a NUB, it is convenient to introduce the 
following reduced volume 

fi 



L n+2 



which is invariant under the overall scaling of the system, L — > j3L and V —> /3 n+2 V. Since 
the volume is conserved and the period L is fixed in our simulations, fi remains a constant 
throughout the evolution. Each class of constant mean curvature surfaces such as that of UBs, 
NUBs, and SDs is parameterized by fi [24]. In particular, the stability analysis of the UB says 
that a UB parameterized by /i is unstable if fi < fi CI = Vt n (^/n/2n) n+1 . 

In the 4D case, the phase diagram of constant mean curvature surfaces [24] tells us that the 
SD, which would be the endpoint of the dynamics, exists even for fi < fi CI with an area smaller 
than that of the UB having the same fi. This suggests the pinch off always happens whenever 
we take fi < fi cr initially, i.e., take 5 > in the initial condition (|T5"|) . We confirmed this by 
taking several values of 5 within the accuracy of our calculation. While in the 12D case, the 
phase diagram [24] tells us that the NUB, which would be the endpoint of the dynamics, exists 
only in 0.28/i cr < fi < fi cr with an area smaller than that of the UB having the same fi. This 
means that if we take a large value of 5 resulting in fi < 0.28/i cr , the perturbed UB cannot 
converge to the NUB but pinches off and approaches a SD, which exists for such a small fi. We 
confirmed this (in the sense described in the previous subsection) by taking several values of 
5 greater than that in Eq. (fT6|) . The parameter e does not play any important role since this 
quantity just changes the dynamical time scale. 

Here, we mention the fate of Rayleigh-Plateau instability in other dimensions within 3 < 
D < 15. For the pair of (e,S) in Eq. (I15p . we observed the pinch off in 3 < D < 11 and the 
convergence to the NUB inl2<D<15, while the non-uniformness of the final equilibrium 
depends on dimension. We should note that if we take 5 too small (0 < S <C 1), the UB 
converges to the NUB even in 10D and 11D since there exists the NUB branch having an area 
smaller than that of the UB [24]. 



S 



4 Discussion 



We have seen that the non-uniform bridge indeed can serve as the global attractor of the 
surface diffusion dynamics associated with the Rayleigh-Plateau instability. Combining this 
result with the correspondence of the phase diagrams between the black hole-black string and 
liquid drop-liquid bridge systems [6,24], it is likely that the fate of the unstable long (but not 
too long) black string is the non-uniform black string around the critical dimension and above. 
On the other hand, although we have seen that the Rayleigh-Plateau instability leads the 
unstable uniform bridge to pinching off in lower dimensions, the pinch off of an event horizon 
will be associated with a naked singularity and other problems peculiar to curved spacetimes. 
Therefore, the relation between the respective endpoints needs further considerations in such 
lower dimensional case^E 

The existence of critical dimensions seems to stem from purely geometrical or morphological 
effects not peculiar to the black strings and minimal surfaces. Therefore, the convergence to 
non-uniform configurations would be realized in other physical systems in which an attractive 
force governs the dynamics. It is interesting to find the convergence to non-uniform configura- 
tions, for instance, in the Rayleigh-Plateau instability in a hydrodynamical system, the Dyson- 
Chandrasekhar- Fermi instability as a Newtonian gravity counterpart of the Gregory-Laflamme 
instability [21,22], and a system of gravitating shell/domain wall. 

One who is familiar with the black hole thermodynamics (e.g., [32]) might wonder that the 
"one-way nature" of the surface area (A < 0) brings to mind the second law of the black hole 
thermodynamics, i.e., the area of an event horizon increases (5Aeu > 0) or the free energy of 
a black hole decreases (5F < 0). Furthermore, that the diffusion dynamics converges to an 
equilibrium state if and only if the mean curvature k becomes a global constant is quite similar 
to the zeroth law, which states the global constancy of the surface gravity on event horizons for 
stationary black holea^). Thus, the surface diffusion dynamics considered in this paper not only 
displays the Rayleigh-Plateau instability similar to the Gregory-Laflamme instability but also 
has some intrinsic aspects each of which has a counterpart in the black hole dynamics. Of course, 
it could not be expected in general that black branes behave like the metal surfaces. However, 
the similarity of phase structures and the correspondences of dynamics will be worth further 
consideration. It is interesting to investigate these correspondences from a more fundamental 
perspective (cf. Ref. [18]). 

As mentioned in Introduction, the behavior near the breakup and a coalescent process in 
the surface diffusion have been known to be described by self-similar solutions [28,33]. The 
self-similar solutions play important roles also in the gravitational collapse, especially in the 
critical phenomena [34,35] and naked singularity formation [36,37]. We saw that depending 
on the parameter p in the initial condition (fT5l) . the unstable uniform bridge will either pinch 
off or converge to the non-uniform bridge. This is analogous to the critical phenomena in the 

2 If the evolution results in pinching off also in the black string dynamics, a certain class of multi-black hole 
configurations should be interpreted as the endpoint of the Gregory-Laflamme instability [31]. 

3 In the theory of surface diffusion, the mean curvature k originates from the chemical potential on metal 
surfaces rather than the temperature [26] , while the surface gravity of a black hole is proportional to the black 
hole temperature. 
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gravitational collapse [34]. It is interesting to investigate the present system from this point of 
view and to look for a counterpart in the black hole-black string system. 



Acknowledgments 

I would like to thank Oscar Dias, Roberto Emparan, Miyuki Koiso, Barak Kol, Kei-ichi Maeda, 
Mu-In Park, and Jennie Traschen for valuable comments and discussions. This work was 
supported by the Golda Meir Fellowship, Israel Science Foundation Grant No. 607/05, and DIP 
Grant H.52. 

A Calculation Details 

We derive the geometric quantities and some identities used in Sec. El 



A.l Geometric Quantities 

A line element of (n + 2)-dimensional flat space in cylindrical coordinates is given by 

g iiv dx il dx v = dR 2 + dz 2 + R 2 d{l 2 n , (19) 

where /i, v — 1, 2, . . . , n + 2 and dQ 2 := a i jdy t dy : > , (i, j — 1, 2, . . . , n) is the line element of a unit 
n-sphere. Let us consider an axisymmetric hypersurface E given by $(-R, z) := R — r(z) = 0. 
Since dR = r'dz on E, we have the induced metric h a b, (a,b = 1, 2, . . . , n + 1) on E and its 
determinant 

h ab dx a dx b = (1 + r' 2 )dz 2 + r 2 dQ 2 n , Vh = r n ^(l + r ,2 )a , (20) 
where a := deta^. The Laplacian on E is given by a formula 

A s = -±=d a (Vh h ab d h ) . (21) 

If we assume that the functions to be differentiated depends only on z, we have A s in Eq. ([S]). 
An outward unit normal of E, denoting it by n^, is given by 

^^$p^ (s "- r ' s;) ' (22) 

The mean curvature can be defined as the trace of an extrinsic curvature of E, which we denote 
by K$. That is, 



k := K» = V/ = gTnu) 



n „ / r' 



R=r(z) rVTTT^ \VT 



r /2 



-Tr=t ■ (23) 



Finally, the normal velocity of E, u, is the projection of the velocity of the surface, which 
we denote by v M , on the unit normal n M , i.e., u := v^n^. Noting that = dtr5 R , we have the 
expression in Eq. ([5]). 
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A. 2 Area Decreasing 

A few identities which might be useful to derive Eq. (jSJ) are presented. Defining 

J : = -7T^ ' (24) 

equation is rewritten as 

r n d t r = -Bd z (r n j) . (25) 

One can easily show an identity, 



dz {vTT^J " v / TT^ + (i + r'2)3/2 • ^ 

Eliminating the second derivative, r", from Eqs. (j4j) and (|26|) . we have 



nr 



n-X. 



n f 



VT+^ = r n K + d z - . (27) 



VT + r' 2 , 

One can show Eq. ([7]) using Eqs. ( l25l) and (127]) via integration by parts. 

B Second-Order Perturbation 

One can extend the linear perturbation in Sec. 12.31 to arbitrarily higher-order ones. First, let 
us expand r(t, z) around the UB, 

oo 

r(t,z) = J2^ r m(t,z) , (28) 

m=0 

where r and ri(t, z) are given in Sec. 12.31 Hereafter, we set r = B = 1 for simplicity. As we 
are concerned only with the backreaction of the linear perturbation ri(t, z), we set r m (0, z) = 0, 
(m > 2). Here, we consider the second-order perturbation. Substituting the expansion ( |28l ) 
into the basic equation, we have the following at 0(e 2 ), 



dtT2 + nr'o + rl 4 ' = n 



\n - 2)r'? + {2ri - r'-l)r'-l - 2r' l rf ) ] . (29) 



After some calculations using the Fourier and Laplace transformations with respect to z and t, 
respectively, we have 



N N 

M (1 - e 2 ^) +J2 B i ( eMkl)t ~ eUl{2kl)t ) cos(2fc^) 

i=l i=l 

+ Yl C ij (e [uJ ^ )+uj{k ^ ]t - e ^ ki+k ^) cos(ki + kj)z 

l<j<j<N 

+ Yl C iJ (e luj{k ' )+uj{k ^ ]t - e^*-"'*) cos(fcj - k 3 )z , (30) 

l<j<j<N 
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where 



_ n 2 R n ( n - 4 - %kf)k1 2 
1 : " ' 4 : " 2[2u{k i ) - u^h)}^ ' 

_ n[(w - 2 t fcjfcjQfcjfcj =f (1 ± fcjfcj)(fc? + k])] 
ij wfc) + ufa) - u{h ± kj) ^ ' ld j 

From the above result, one recognizes an interesting non-linear effect. If we take fc, > \fn 1 
(1 < i < iV, note that we set ro = 1), the linear perturbation and almost all terms in the 
second-order perturbation (I3UI) are suppressed exponentially or remain to be constant at most. 
However, the term containing factor e ul ( ki ~ kj ' )t in Eq. (15U1) grows if there is a couple of sinusoidal 
waves satisfying \ki — kj\ < y/n, (i ^ j) because uj{ki — kj) > for such a couple. This fact 
leads to the interesting behavior that high wavenumber modes begin to grow after an initial 
decay followed by an incubation time [27]. 
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